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We study the density of states and localization properties of the lowest Landau levels of graphene 
at high magnetic fields. We focus on the effects caused by correlated long-range hopping disorder, 
which, in exfoliated graphene, is induced by static ripples. We find that the broadening of the lowest 
Landau level shrinks exponentially with increasing disorder correlation length. At the same time, 
the broadening grows linearly with magnetic field and with disorder amplitudes. The lowest Landau 
level peak shows a robust splitting, whose origin we identify as the breaking of the sublattice (valley) 
degeneracy. 

PACS numbers: 73.43.-f,73.63.-b,71.70.Di 



I. INTRODUCTION 

The observation of the anomalous quantum Hall 
effect 1 2 is one of the most striking and robust mani- 
festations of the underlying massless Dirac fermions in 
graphene near half filling. The energy scales in graphene 
are such that the quantization of the Hall plateaus can be 
observed even at room temperature at sufficiently high 
magnetic fields. 3 The energetics in graphene also favors 
a direct experimental access to the low-lying Landau lev- 
els by infrared spectroscopy^!! and by scanning tunnel- 
ing spectroscopy^^ something hardly possible in conven- 
tional semiconductors. 

Further information about the nature of the lowest 
Landau levels in graphene has been recently obtained 
in thermal activation experiments ptMl showing that the 
zeroth Landau level (n = 0) is much sharper than the 
first and higher Landau levels. These observations are 
the main motivation of the theoretical study presented 
in this paper. 

Disorder is key to understand the electronic transport 
properties in graphenepK^ particularly in the quantum 
Hall regime, where the conductivity plateaus are con- 
ventionally explained by delocalized states surrounded 
by localized ones. However, the mechanisms that lead 
to localization in QHE in graphene are still not clear.^ 
Currently, there is still some debate on the most rele- 
vant disorder mechanisms for transport in graphene.^ 
Among those, ripple disorder is believed to play an im- 
portant role. Static ripples give rise to random correlated 
hopping disorder] 19 * 20 * which is the disorder mechanism 
analyzed in this paper. 

The shape and width of the lowest Landau levels (LLs) 
in graphene have been investigated in several theoretical 
studiesPf^ The broadening of the LLs differs among 
disorder models. In particular, the inclusion of a finite 
correlation length on the hopping (off-diagonal) disorder 
model was recently reported to induce an anomalously 
sharp n = LL compared to higher levels PUSH it wa s 



found that the width of the zeroth LL (To) shrinks to 
zero as soon as the hopping correlation length A exceeds 
the lattice parameter a, in line with analytical studies of 
the effects of long-range chiral disorder.^ 

Regarding the localization properties of the low- 
est LL, numerical simulations using uncorrelated hop- 
ping di sordeJ ^ 2 ^ and white- noise random magnetic flux 
disordeiJ^ 3 -^ observe an interesting distinct qualitative 
feature in the quantum Hall spectrum of graphene, 
namely, a splitting. It was found that, in such chiral 
disorder models, the lowest LL splits into two Gaussian 
shaped peaks, even in the absence of both a Zeeman term 
and electron-electron interactions. The splitting energy 
AE is linearly proportional to the disorder strength and 
scales with th e squa re root of the applied perpendicular 
magnetic field! 24 * 31 ! A similar square root magnetic field 
dependence of the splitting of the n = Landau level has 
been experimentally observed in Ref. |9l 

In this paper, instead of the white-noise random hop- 
ping (or magnetic flux) model, we address the more real- 
istic correlated random hopping disorder modelP^^ We 
present a systematic study of the shape of the lowest 
LLs and their localization properties as a function of the 
hopping disorder correlation length, as well as of other 
relevant parameters of the system, such as disorder am- 
plitude and magnetic field. We find that Tq decays expo- 
nentially with the correlation length, never fully vanish- 
ing for any finite A. More importantly, we observe that 
the ratio Ti/Tq depends only on the disorder correlation 
length, showing no significant variation neither with the 
disorder strength, nor with the magnitude of the applied 
magnetic field B, provided the system is in the quantum 
Hall regime. In addition, we study the splitting of the 
n = LL, which is inferred from the analysis of the par- 
ticipation ratio. We show that this splitting shrinks with 
increasing values of A, but is still present even for cor- 
relation lengths for which the n = LL width becomes 
very small. 

The paper is organized as follows. In Sec. [Hj we 
present the model used in our numerical simulations. 
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The analytical framework for the interpretation of our b ^ 34 * 35 * 



results is discussed in Sec.|III| Next, we analyze the spec- 
tral, Sec. |IV[ and localization properties, Sec. [V| of the 
model. We conclude summarizing our results and dis- 
cussing their relevance to the interpretation of experi- 



ments on the quantum Hall effect in graphene in Sec. VI 



II. MODEL DESCRIPTION 

Graphene is a monolayer honeycomb lattice of carbon 
atoms with a lattice constant a = 2.46 A. Its primitive 
unit cell contains a pair of atoms that form two triangu- 
lar sublattices, denoted by A and B. The tight-binding 
Hamiltonian model for a graphene monolayer reads 
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where the sum runs over nearest-neighbor sites. The ex- 
ternal magnetic field £?, perpendicular to the graphene 
sheet, is included by Peierls' substitution, namely, <pij = 
2ir(e/h) J^ z A-dl . In the Landau gauge, A = (0,Bx) 
and considering a brick wall lattice, which is topologi- 
cal^ equivalent to the hexagonal lattice,^ one has <pij =0 
along the x direction and faj = ±7r(x/a)$>/$>o along the 
=F?/ direction, with $/<^> = Ba 2 ^f3e/(2h) per unit cell. 

The random hopping disorder is implemented by ran- 
domly choosing the hopping parameters Uj from a uni- 
form distribution of width W around the average value 
t = 2.7 eV. In addition, we impose here a spatial cor- 
relation to the hoppings following a Gaussian profile of 
width A. Figure [I] illustrates typical realizations of the 
disordered hopping parameter tij for two different corre- 
lation lengths (A = 2a, 4a). The color scale refers to the 
hopping amplitude at the middle point between the two 
nearest-neighbor sites i and j. As expected, Fig. [I] shows 
smoother (less abrupt) variations of the hopping values 
with increasing correlation length. 

We consider graphene lattices of MxN carbon atoms 
(M zigzag chains, each containing N atoms) with peri- 
odic boundary conditions. The linear sizes of the lattices 
are 
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Most of the numerical results shown in this paper are 
calculated for lattices of MxN = 100x90 atoms, corre- 
sponding to a size L x xL y = 21.1 nm x 10.9 nm, the 
same lattice dimensions shown in Fig. [I] 



III. DIRAC HAMILTONIAN: ANALYTICAL 
RESULTS 

Near half filling, the low-energy properties of the tight- 
binding Hamiltonian in Eq. ([I]) are described by noninter- 
acting massless Dirac fermions in an uniform perpendic- 
ular magnetic field, with an effective Hamiltonian given 
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(3) 

where Hvp = V3at/2 is the Fermi velocity, 7r = p+eA/c, 
and p stands for the electron momentum operator. The 
Hamiltonian ([3| operates on a four-components wave 

function (^a^b^a^b')i wnere and ®b re P" 
resent the envelope functions at A and B sites for the K 
point and and Sfrg for K 1 ', respectively. In this pa- 
per we do not consider explicitly the electron spin degree 
of freedom and all states are assumed spin degenerate. 

The eigenstates of the Hamiltonian (|3| are labeled by 
(j, n, fc), with the valley index j = K^K', the Landau 
level index n = 0, ±1, ... , and the wave vector k along 
y direction.^ The eigenenergy depends solely on n as 
e n = hbJB sgn(n)Y / 1^J, where Hub = V2j/£b with the 
magnetic length given by £b = yh/eB. For ^/^o < 
0.05, which implies £b/cl ^> 1, lattice size effects have a 
negligible influence on the graphene LLs in this model. 
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FIG. 1: Typical spatial "landscapes" of the fluctuation in 
the hopping matrix elements SUj /W = (Uj — t)/W for two 
different correlation lengths (A = 2a, 4a). The color scale rep- 
resents the magnitude of 8Uj/W. The lattices shown above 
have MxN — 100x90 atoms, which corresponds to lateral 
dimensions L x xL y — 21.1 nm x 10.9 nm. 
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The eigenfunctions are written as 
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Here C n = 1 for n = 0, C n = 1/V2 for n^O, and 
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with z = (x — kl 2 B )/lB and H n (z) denoting Hermite poly- 
nomials. 

It was realized early 2 ^ that the level with n = is 
special since its amplitude is non-zero only in one of the 
sublattices, namely, at B sites for K and A sites for K' . 
Consequently, while a random on-site disorder potential 
gives only intravalley mixing within either the K and K' 
valleys, random hopping causes intervalley mixing. (No- 
tice that this is quite the opposite of what occurs at zero 
magnetic field when diagonal disorder is present.) The 
wave function in LLs with has nonzero amplitudes 

on both A and B sites, so that intervalley mixing is al- 
ways possible. 

In this study, we consider hopping disorder caused by 
randomness in the hopping integral connecting neighbor- 
ing A and B sites. This disorder can be long ranged, 
when caused by static ripples, or short ranged, when orig- 
inated by scatterers located at points in-between neigh- 
boring sites. 

Assuming that t shifts to t + St between neighboring 
sites Ra and hopping disorder gives rise to a short- 
range potential given by^ 
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with u h (V3a 2 /2)5t, z x = e* KR *, and z' x 
e %YL n x f or x = A and B. For a hopping disorder 
concentration rih, the self-consistent Born approxima- 
tion estimates the Landau level broadening as r s h or t = 
(2nh) 1 ^ 2 Wh/( / 7T^ J B), independent of Landau level index. 
Numerical simulations 24 exhibit the same scaling of 
r s hort with V^5, but show that r s h or t increases with the 
index n. 

Alternatively, when the randomness in the hopping in- 
tegrals shows long-range correlations, the disorder Hamil- 
tonian can be form ulated in terms of an effective random 
magnetic field! 2Q | 36 ^ In this case, one assumes that lattice 
deformations cause a smooth shift in the hopping inte- 
grals between the site j and its three nearest neighbors i. 
At low energies, this effect can be incorporated into Dirac 



Hamiltonian by introducing an effective vector potential 
that reads 



(8) 



Here, Af = Af =F iAf, e 7 are vector connecting 
a lattice site to its neighbors, and q w ±qr , where 
qo = 47r/(3v // 3«)e 2/ with the Cartesian coordinate x run- 
ning along the armchair direction. The subscript + (— ) 
corresponds to the K (K') valley. Random hopping is 
accounted for by adding the effective vector potential 
A eff to the momentum operator 7r appearing in Eq. (|3|, 
namely} 3 -^ 
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Notice that the effect of A*± is to locally shift the Dirac 
cones K and K' in opposite directions and there is no val- 
ley mixing. The structure of H± immediately reveals that 
the n = states are unique: Since (^^ k \H — Ho\^/^ k ,} = 
0, in lowest order, long-range hopping disorder does not 
affect the n = states. 

Unfortunately, there is no such simple picture for 
treating the crossover regime between short- and long- 
range hopping disorder. In the following, we interpret 
the results of our numerical simulations by invoking the 
pure long-range description provided by Eq. ([9|, what is 
known about short-range hopping disorder, and by build- 
ing a plausible interpolation between these two limits. 



IV. SPECTRAL PROPERTIES 

In this Section, we analyze the spectral properties of 
graphene obtained from the tight-binding model with 
correlated random hopping presented in Sec. [TTJ We 
will focus our attention on the width of the disorder- 
broadened Landau levels and its dependence on the dis- 
order correlation length A, on the disorder amplitude VF, 
on the magnetic flux <E>, and on the LL index n. 

Figure [2^ shows the density of states (DOS) corre- 
sponding to the four lowest LLs (n = 0,1,2, and 3) 
broadened by a correlated random hopping disorder of 
amplitude W/t = 0.3 for different correlation lengths, 
namely, A/a = 1, 2, and 4. Since particle-hole symme- 
try is preserved by the nearest-neighbor hopping disorder 
model, we only show the n > LLs. The magnetic flux 
considered is $/$ =0.02. The results typically corre- 
spond to averages over 600 disorder realizations. 

Figures [2]o and ^ are zooms of the DOS around the 
n=0 and the n=l LLs, respectively, indicating that the 
broadening shrinks for increasing values of A. An inspec- 
tion of Fig. [2J3 shows that, for the n = 1 LL, an increase 
in the correlation length causes only small modifications 
on the shape of this Landau band, plus an overall reduc- 
tion of the width Y\. Higher Landau levels, n > 1, show 
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FIG. 2: Density of states (in arbitrary units) for a magnetic 
flux $/$o=0.02 for three different correlation lengths A/a of 
the random hopping potential (in all cases W/t = 0.3 and the 
lattice size is MxN— 100x90 sites), (a) DOS showing the 
n=0, 1, 2, and 3 Landau bands, (b) Zoom of the n — LL 
states, (c) Zoom of the n — 1 LL states. 



transition to zero width was observed. It is worth men- 
tioning that, when A increases, making the width of the 
n = LL much smaller than the width of the higher lev- 
els, it is important to use DOS histograms with a much 
finer energy resolution around the n = LL than for the 
higher LLs. Not doing so can lead to an erroneous im- 
pression that the n = LL width vanishes sharply as the 
correlation length increases. 

In the following, we describe how we define and quan- 
tify the width T n of the disordered-broadened LLs. For 
the n > LLs, which display a Gaussian-like shape, T n 
is taken as the full width at half- height. Figure [2] clearly 
shows that the DOS shape of the n = LL is very dif- 
ferent from the higher LLs. As pointed out in Ref. [24j 
an off-diagonal disorder model induces a splitting of the 
n = LL into two degeneracy-broken n = Landau 
bands, causing the observed DOS shape for the n = LL, 
namely, not fully split levels. This is so because the en- 
ergy splitting always has the same order of magnitude of 
the LL broadening. Here, the n = LL shapes observed 
in Fig. |2j) can be reasonably well fitted by a superposi- 
tion of two equal Gaussian curves. Therefore, the n = 
LL width is considered as the width at half-height of one 
of these superposed bands (which is approximately half 
the width at half- height of the entire band). 

Figure [3] shows how the LL widths decrease as the hop- 
ping disorder correlation length grows. In Fig. [3^l, one 
observes that not only r decreases with A, but all other 
LLs do so. When examining the same data in a log- 
linear graph (Fig. |3|d), one observes that the n = LL 
behaves quite differently from the others. For n > 0, the 
widths T n decrease slowly with increasing A/a. Figure 
[3] shows an apparent tendency to saturation at a value 
T n /t « 0.01. We computed r n> o for larger values of A/a 
(not shown here) and concluded that this in not quite cor- 
rect. Instead, we observe that the rate by which all r n>0 
decrease becomes smaller as A /a becomes larger. The be- 
havior is very different for the n=0 LL, whose width de- 
cays exponentially with A. This decay is sustained down 
to the numerical precision of our simulations. 

Further insight about the effect of long-range corre- 
lated hopping disorder on the levels widths Y n can be 
gained from a perturbation theory analysis. Let us de- 
note the disorder gauge potential by V = H± —Hq, where 
both H and H are defined in Sec. |III| 

In first order, the matrix elements required to calcu- 
late the energy corrections for the degenerate states that 
belong to the nth Landau level at the K valley are 



the same qualitative behavior as n = 1. In contrast, for 
the n=0 LL (Fig. [2}d), one observes a much stronger 
suppression of the level broadening upon increasing A, 
in agreement with numerical results obtained in Ref. l25l 
However, the inset of Fig. (2}d shows that, despite the 
large reduction of To with increasing A, the width never 
goes to zero within the range of A we used. This is at 
odds with the results reported in Ref.[25l where an abrupt 
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Since long-range hopping disorder does not mix valleys, 
K is a good quantum number in this model. As dis- 



cussed in Sec 



qua 
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0. For n > 0, the situation is 



different. Exact diagonalization in a (n > 0, K) subspace 
involves matrix elements of the kind V^ K . kk >- Due to 
the spinor structure of the evaluation of such ma- 
trix elements amounts to the spatial integration of the 
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FIG. 3: Width of the nth Landau level T n /t, n = 0, 1, 2, 3, 4, 
as a function of the disorder correlation length X/a for $/$o = 
0.02 and W/t = 0.3. (a) Linear scale, (b) Log- linear scale. 



product exp[z(fc — k')y)(j) n k> {x)(j) n -i^{x) times the Gaus- 
sian fluctuating gauge potential. This results in non-zero 
matrix elements, but they are very quickly suppressed as 
X/£b becomes large. 
Let us define 
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to help us to discuss second-order effects. The evaluation 

°f Vn=o K-kk' i nv °l ves matrix elements where products of 
and 0o, k appear. For long-range disorder, such ma- 
trix elements vanish with increasing X/£b- In summary, 
the perturbation fails to mix states within the n = 
multiplet and also with n' ^ 0. Hence, we expect the 
same behavior at all order of perturbation theory. This 
is consistent with statement that Tq = for long-range 
hopping disorder P^Using the reasoning presented above, 
matrix elements of the kind V^J Q K . kk , involve, among 
other components, the integration of products wave func- 
tion amplitudes such as exp[z(fc — k')y](j) n k'(x)(j) n ^k{x). 
These become small for k ^ k' in the limit of X/£b ^> 1 
and, for k = k' , are quite independent of X/£b ^> 1- 

This analysis rules out long-range hopping disorder as 
the mechanism behind the exponential suppression of Tq 



with increasing A. We speculate that this behavior is 
caused by matrix elements that admix valleys, a remnant 
of the crossover regime. As a consequence, we expect Tq 
to scale linearly with the disorder strength W and the 
eigenstates to be a superposition of a wave functions with 
probability amplitudes in both sublattices. Our simula- 
tions are in line with the latter statement, as we discuss 
below. 

When a magnetic field is present, there is an impor- 
tant length scale to be considered, namely, the mag- 
netic length; for convenience, let us write it in the form 
£b/cl=0.371/ \/&/&o. For the magnetic flux used in our 
simulations (<£/<I>o=0.02), we obtain £b / a=2.62, which 
is close to the values of X/a used as well. Therefore, we 
expect the magnetic length to play a role in any inter- 
pretation of our results. Indeed, Fig. [3^i suggests that 
the LL widths change its dependence with A for A > £b- 
However, Fig. (3}d makes clear that a slow dependence on 
A only occurs for the n >0 LL. This is in line with our 
perturbative analysis. Unfortunately, the picture is not 
entirely consistent: We expect the second-order terms to 
be dominant in the calculation of the broadening of r n> o 
for X/£b : which is not observed in the simulations (see 
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FIG. 4: (a) Landau level width T n /t as a function of the hop- 
ping disorder strength W/t for the n=0 and n=l LLs (mag- 
netic flux fixed at <£/<i>o=0.02 and A/a=1.5). (b) F n as a 
function of magnetic flux $/$o for the same LLs (now for a 
fixed disorder strength W/£=0.3 and A/a=2.0). 
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FIG. 5: (a) Ratio between different Landau level widths 
r n /r n _i as a function of A/a. The ratio between higher LLs 
widths, like Y2/Y1, is quite independent of the correlation 
length A. In contrast, the ratio Y\/Yq grows rapidly with A/a 
(notice the logarithmic scale), (b) Ratio Yi/Yo as a function 
of the disorder strength W/t for A/a = 1.5. (c) Ratio Y\/Yo 
as a function of the magnetic flux <f>/<I>o for A/a = 2.0. The 
panels (b) and (c) show that Yi/Yq is almost independent of 
W/t and $/$o- In the log-linear plot of panel (a) the results 
corresponding to different values of W/t and $/$o collapse 
under the same dot. 



discussion below). This remains to be understood. 

We call attention to the fact that the authors of Ref. [25] 
used the same flux value we considered here, as well as 
a similar range for the correlation length A. However, in 
our calculations, a finite width of the n = LL can be 
seen even for A/a=4 within the numerical precision we 
use. We checked (not shown here) that these results are 
not influenced by varying the system sizes, i.e, there are 
no finite lattice-size effects in the parameter region we 
investigated. 

We have also investigated the dependence of the n = 
and n = 1 LLs widths on the disorder and magnetic field 
amplitudes (Fig. [4|. For both parameters, there is a 
clear linear increase. 

Since the LL widths depend linearly on W and on 
$/$o, we conclude that there is universality in the behav- 
ior of the Landau level widths, namely, the ratio between 
different Landau level widths, Y n /Y n ^ depends solely on 
A/a. This is illustrated in Fig. [5^i, where one can see the 
ratio Yi/Yq growing rapidly with A/a (notice the loga- 
rithmic scale), while Y2/T1 remains essentially constant. 
This novel result allows our simulations to make contact 
with the experiments. Notice that due to computational 
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FIG. 6: Landau level width Y n /t as a function of LL index n 
for three different values of A/a, considering $/$o = 0.02 and 
W/t = 0.3 (same magnetic flux, disorder amplitude, and lat- 
tice sizes considered in Fig. |3j). The dashed line corresponds 
to a fitting with square root dependence, which fits well only 
the A/a = 2 case (see text). 



limitations, our lattices sizes constrain us to consider un- 
realistically large values of magnetic field magnitudes. 
However, since the ratios of LL widths are rather insen- 
sitive to the values of W and ^/^ch as illustrated by 
Figs. [5)3 and c, we expect this result to apply for realistic 
settings as well. 

The dependence of the width on the LL index is shown 
in Fig. [6j In this case as well, Y n increases with n when 
A, W, and are fixed (notice that the values for these 
parameters are the same considered in Fig. [3|. We find 
that for A/a = 2, the Y n versus n curve is perfectly fitted 
by the functional form y = A^fx, with A=0.0107 (dashed 
line in Fig. |6|. For other values of A, including the cases 
A/a = 1 and A/a = 4, the numerical data do not show a 
square root dependence on LL index n. At this point, it is 
worth comparing these results to the effect of a diagonal 
disorder on the DOS in graphene at the quantum Hall 
regime. The Landau level broadening Y n for a diagonal 
w hite-n oise disorder is quite independent of the LL index 
n p 2 * 23 1 while for the correlated diagonal disorder, Y n is 
observed to slightly decrease with increasing n. 22 These 
effects for the diagonal disorder in graphene are similar 
to the observed for conventional quantum Hall systems 
with diagonal disorder modelsp^but are in clear contrast 
to the increase of Yr, with n observed here. 



LOCALIZATION PROPERTIES 



While in Sec. [IV] we considered the DOS and analyzed 
the LL widths as a function of several parameters, we now 
investigate the localization properties of states within the 
LLs. To infer the degree of localization of the states we 
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use the participation ratio (PR), which is defined as^ 



PR 



(12) 



where ipi is the amplitude of the normalized wave func- 
tion on site i and N' = M x TV is the total number of 
lattice sites. The PR is therefore directly related to the 
proportion of the lattice sites over which the wave func- 
tion is spread: the PR for a localized state vanishes in 
the thermodynamic limit, while peaks in the PR indicate 
the presence of extended states (critical energies). 

Figure [7^ shows the calculated PR in an energy win- 
dow comprising the Landau levels n=0, 1, 2, and 3. 
While the PR for the n > LLs indicates the presence of 
localized states at Landau band tails and extended states 
at the band middle, as expected, the PR of the states 
from the n = LL shows a double hump structure, i.e., a 
splitting into two peaks. The splitting is more clearly ob- 
served in the energy scale zooms of Fig. [TJd and c. This 
splitting of critical energies in the lowest LL was pre- 
viously observed and discussed for uncorrelated random 
hopping disorde r 21 ! 24 * (and for the sim ilar c ase of uncor- 
related random magnetic flux disorder^^ 2 -') . The novel 
and interesting aspect we observe here is that the split- 
ting is rather robust and survives even the sharp width 
reduction of the n=0 LL due to the increasing correla- 
tion length. While previ ous w orks considered correlated 
random hopping disorder p£126| they did not calculate the 
localization properties of the states within LLs and there- 
fore missed this splitting of critical states at n = 0. 

Although experiments 9 have observed this splitting en- 
ergy in the n=0 LL, it is hard to make a quantitative 
comparison due to the lattice size compromises we need 
to make in our simulations. Moreover, we believe that a 
full explanation of the experimentally observed splitting 
requires taking into account the Zeeman term and possi- 
bly electron-electron interactions, which calls for further 
theoretical investigation. 

Due to particle-hole symmetry, the probability densi- 
ties \^(E)\ 2 and the PR of states at energies — E and 
+E are identical (see, for instance, Ref. [16| and refer- 
ences therein). However, when looking directly at the 
wave function amplitudes ^(E), one can see a differ- 
ence between the split states at — E and +E: In one 
of the sublattices (sublattice A, for example) the am- 
plitudes are exactly the same, while in the other sub- 
lattice (B) they have the same magnitude but oppo- 
site signs, that is, ^ A (~E) = "&a(+E) and V B (-E) = 
— ^b(+E)- This guarantees that conjugated particle- 
hole states are orthogonal and have the same probability 
densities \^(E)\ 2 . It is noteworthy that orthogonality 
imposes that the probability weight at both sublattices 
is the same. The split states are therefore similar to 
bonding-antibonding states. The table in Fig. [8] summa- 
rizes these features. 

Another feature observed in Fig. [7] is the effect of in- 
creasing correlation length. The effect on the n=0 LL 
is completely different (opposite) to what is observed in 
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FIG. 7: Participation Ratio as a function of E/t for different 
random hopping correlation lengths A/a. Here ^/^o = 0.02, 
W/t = 0.3, and MxiV=100x90 (namely, the same magnetic 
flux, disorder amplitude, and lattice sizes considered in Fig. 
|3j). (a) Energy window containing the four lowest LLs, n—{). 1, 
2, and 3. (b) Zoom on the n — Landau level, (c) Amplified 
zoom on the n — states. 



higher levels. In higher levels, the increase in A causes an 
overall reduction of PR values, while in the n=0 LL, de- 



spite the suppression of broadening (see Sec. IV), one ob- 
serves an overall increase in the PR. This difference in be- 
havior is probably related to the fact that the states from 
the higher levels have much longer localization lengths 
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FIG. 8: Properties of the wave functions ^ for n — LL states 
at energies +E and — E. While the PR and the probability 
densities of the wave functions |^(^)| 2 are identical for states 
at energies — E and +E, the wave function amplitudes ty(E) 
show differences in only one of the sublattices, in analogy to 
bonding- antibonding states (see text). 



than those from the states with n=0. 3T To investigate 
possible finite-size effects on the n = LL when A/a 
is large, we ran simulations with larger lattice sizes (not 
shown). We found no change in the position of the identi- 
fied critical states, as well as no change in the LL widths. 



VI. CONCLUSIONS 

We have investigated the effects of spatially correlated 
random hopping disorder on the structure of Landau lev- 
els in graphene. We quantified the behavior of the nth 
Landau level width T n as a function of the correlation 
length, as well as other relevant parameters of the system: 
disorder amplitude, magnetic field, and Landau level in- 
dex n. We found that T n gets narrower with increasing 
correlation length for all Landau levels, and not only for 
the n = level. However, a logarithmic plot of T n as 
function of A (Fig. 3b) clearly showed that while for 
n > the widths decrease slowly, for n = they decay 
exponentially with increasing A. We found no sign of any 
abrupt vanishing of the n = Landau level width at finite 
A. This suggests that a different physical mechanism is 
behind the narrowing of the n = LL when compared to 
all the other levels (specially when the correlation length 
A becomes of the same order of magnitude or higher than 
the magnetic length £b)- We speculate that the under- 
lying mechanism is due to valley mixing, remnant of the 
crossover regime. 



T n increases linearly with both magnetic field and dis- 
order amplitude. Another interesting observation is that 
for any fixed correlation length, T n always increases with 
increasing LL index n, which is completely different to 
the dependence observed for diagonal (on-site) disorder 
models. More importantly, we observe that the ratio be- 
tween the n=l and n=0 Landau level widths, Ti/Tq^ de- 
pends only on the correlation length A and is rather in- 
sensitive to the disorder strength and to the magnitude 
of the applied magnetic field. This allows a closer con- 
tact of our results with experiments. In Ref. [TTJ for 
instance, the authors experimentally observe a n = OLL 
width of about 20K, while the width of higher levels is 
observed to be about 400K that is Ti/Tq = 20. Using 
this information and the results of our Fig. [5j we infer 
that A/a = 2.2. This estimate gives a lower bound for 
the corrugation length of graphene in the experiment of 
Ref. |TTJ since it neglects all disorder broadening sources 
but hopping disorder. 

We also considered the role played by changing the 
hopping disorder correlation length on the localization 
properties of the states within the different Landau levels. 
The splitting of two critical energies for the n = level, 
previously reported for uncorrelated random hoppingsp^ 
is still clearly defined for correlated hopping disorder. It 
is a robust effect even with the sharp width reduction of 
the n — level that occurs for large values of A. The 
n = level can therefore be considered as a superposi- 
tion of two levels with broken degeneracy, symmetrically 
split around E = 0. In addition, after analyzing the wave 
functions of states belonging to the n = level, we were 
able to identify a symmetric structure of these states. We 
found that although any two states at — E and +E have 
the same probability density and also the same participa- 
tion ratio, in one sublattice the wave function amplitudes 
of states with energies ±E are exactly the same, while 
in the other sublattice these states have the same mag- 
nitude but opposite signs. Therefore, the origin for this 
n = level splitting is clearly related to the breaking 
of the sublattice degeneracy induced by the hopping dis- 
order, a further manifestation of the influence of valley 
mixing disorder. 
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